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In this paper we study the renormalization of the Schwinger-Dyson equations that arise in the 
auxiliary field formulation of the O(N) (j> 4 field theory. The auxiliary field formulation allows a 
simple interpretation of the large-N expansion as a loop expansion of the generating functional in the 
auxiliary field \, once the effective action is obtained by integrating over the <f> fields. Our all orders 
result is then used to obtain finite renormalized Schwinger-Dyson equations based on truncation 
expansions which utilize the two-particle irreducible (2-PI) generating function formalism. We first 
do an all orders renormalization of the two- and three-point function equations in the vacuum 
sector. This result is then used to obtain explicitly finite and renormalization constant independent 
self-consistent S-D equations valid to order 1/N, in both 2+1 and 3+1 dimensions. We compare 
the results for the real and imaginary parts of the renormalized Green's functions with the related 
sunset approximation to the 2-PI equations discussed by Van Hees and Knoll, and comment on the 
importance of the Landau pole effect. 
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I. INTRODUCTION 



Recently there has been interest in studying field the- 
ory using two-particle irreducible (2-PI) methods Q in 
both finite temperature IS E1FJ and non-equilibrium sit- 
uations [1 H Q, 11 @, El Hi IH III G3 The value 
of the 2-PI formalism for non-equilibrium problems is 
that it allows one to make approximations that go be- 
yond the Hartree or large-N approximation without en- 
countering the serious problems of secularity found in a 
straightforward expansion about the Hartree or leading- 
order large-N approximation using the generating func- 
tional or. equivalently, the one-particle irreducible (1-PI) 
action The 2-PI methods lead to self-consistent 

equations for the Green's functions which require non- 
perturbative renormalization. Recently the renormaliza- 
tion of the equations obtained from the 2-PI approach 
applied to the standard formulation of (f) A field theory 
(first discussed by Calzetta and Hu [ll|) has been con- 
sidered by both Van Hees and Knoll LjOf . and by Blaizot 
et. al. . This direct loop expansion of cj) 4 field theory 
is a summation of the coupling constant expansion and 
needs to be resummed in order to be related to a 1/N 
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expansion The approach to renormalization in the 
above works was based on formal Bogoliubov-Parasiuk- 
Hepp-Zimmermann (BPHZ) [3 and dimensional regu- 
larization Il9| methods, rather than multiplicative renor- 
malization |2l |^| . The advantage of the multiplica- 
tive renormalization approach for initial-value problems 
is that it lends itself more easily to the momentum space 
cutoffs that occur when one uses numerical methods to 
solve the integro-differential equations of the closed time 
path (CTP) formalism |22|. In non-equilibrium situa- 
tions, the Green's function equations are usually only 
spatially translational invariant and to make the calcula- 
tional tractable, a maximum 3-momentum is introduced 
(3-momentum cut-off A). In dynamical situations the 
calculational schemes that have been used usually rely 
on mode expansions for the quantum fields which intro- 
duce non-covariant momentum cutoffs. Therefore, for 
practical reasons it is useful to consider direct renormal- 
ization of the Schwinger-Dyson (S-D) equations that are 
rendered finite by momentum space cutoffs. Such an ap- 
proach was very useful in the leading order in large-N ap- 
proximation, where we found that using lattice versions 
of the renormalization scheme gave us results (as long as 
we were far from the Landau pole) that became indepen- 
dent of the cutoff for a wide range of cutoffs when we 
kept renormalized parameters fixed j^l |^| • It was also 
important when studying the time evolution of the (ex- 
pectation value of the) energy-momentum tensor to un- 
derstand the non-covariant nature of the cutoff scheme, 
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so that the correct physical energy densities and pres- 
sures could be extracted from the (non-covariant) situ- 
ation arising in the truncation scheme used for numeri- 
cal simulation [25| . By not automatically subtracting off 
the logarithmical (log) divergences related to coupling 
constant renormalization, we were able to have another 
check on the numerical simulations by studying how the 
simulations became independent of the cutoff for fixed 
renormalized parameters. 

In any truncation scheme, such as expanding the 2- 
PI generating functional in terms of loops or 1/N, the 
renormalization has to be guided by the structure of the 
exact renormalized S-D equations. Thus as a preliminary 
step to renormalizing the truncated S-D equations, one 
needs to know the structure of the exact renormalized 
S-D equations. The strategy for obtaining the renormal- 
ized S-D equations and then using them to renormalize 
the next to leading order in 1 /N is discussed in [2l], |^| • 
In those papers, however, the perturbative 1/N expan- 
sion was discussed and not the resummcd 1/N expansion 
obtained from the 2-PI formalism. 

Our procedure is as follows: First we derive the exact 
S-D equations for the auxiliary field formalism. These S- 
D equations are simpler than those of the original formu- 
lation of (/) 4 field theory because the only quantities that 
need renormalization are the propagators for the field tfi 
and the auxiliary field x, as we U as the three-particle 
X<W> irreducible vertex. Analyzing graphs one finds that 
once those quantities are renormalized one never gener- 
ates a new divergence in the coupling of four 4> particles. 
In the cutoff S-D approach to renormalization, one has 
to show not only that the renormalized equations are 
finite, but also that they are independent of all the (in- 
finite) renormalization constants. That is, all the equa- 
tions need to be written in terms of the renormalized 
Green's functions, vertices and masses. The reason for 
using the auxiliary field formulation of A</> 4 field theory 
is that the 1/N expansion has a simple interpretation as 
a loop expansion [in the auxiliary field] of the generating 
functional of the effective action obtaine d by integrating 
out the scalar fields keeping x constant |2fl |27| • The S- 
D equations arising from this auxiliary field formulation 
was first discussed in Ref. psj . 

The 1 /N expansion is an asymmetric expansion which 
treats the field exactly (with x fixed), and then counts 
loops in x- Thus the 2-PI formalism, which treats <f> 
and x on an equal footing, is not a natural formalism 
for incorporating the large-N approximation except in 
leading order. Its main virtue is that it leads to self- 
consistent approximations that are energy conserving 
and non-secular when applied in non-equilibrium con- 
texts. The basic propagators that occur in the large- 
N expansion, when viewed as propagators coming from 
the effective Lagrangian obtained after integrating out 
the (j> field in the path integral, have different behav- 
ior with regard to 1/N @: the <p(f> propagator G being 
0(1), the XX propagator D being of order 1/N and the 
x4> mixed propagator K (which vanishes if symmetry is 



unbroken) is also of order 1/N. Thus the first nontriv- 
ial two- loop 2-PI vacuum graph has terms GGD ~ 1 /N 
and KKG ~ 1/N 2 . At the level of the equation for 
the inverse two-point Green's functions, the first non- 
trivial approximation (counting loops) has no vertex cor- 
rections, but mixes orders of 1/N. Two approximations 
which have been studied recently in studies of thcrmal- 
ization have been based on keeping one or both of these 
two- loop graphs 0, Q • The first approximation has been 
called the 2-PI 1/N expansion and the second the bare- 
vertex approximation (BVA). Both these approximations 
are identical when (</>) = and thus the renormaliza- 
tion scheme is identical for both approximations. We 
find that both approximations after renormalization re- 
quire that the renormalized vertex function is set to 1. 
Thus the BVA is a misnomer in dimensions greater than 
2 + 1 where wave function and vertex function rcnormal- 
izations are necessary. 

Recently, the 2-PI 1/N expansion has been used to 
study the nonequilibrium dynamics of field theory in 3+1 
dimensions: First, in Ref. the 2-PI 1/N was used to 
study the parametric resonance of an O(N) symmetric 
scalar theory, at very weak coupling constant. Secondly, 
in Ref [Toj , the same approximation was used to investi- 
gate the nonequilibrium dynamics of a 3+1 dimensional 
theory with Dirac fcrmions coupled to scalars via a chi- 
rally invariant Yukawa interaction. In the later case, the 
system was shown to reach, at late times, a state which 
was well described by a thermal distribution. In the 
above work, only the quadratic mass divergences were 
eliminated, and the renormalized mass extracted from 
the oscillation frequency. Determining the renormalized 
mass dynamically makes it a bit more difficult to study 
the behavior of the theory for fixed renormalized mass as 
a function of the momentum space cutoff. The advantage 
of having renormalized S-D equations which depend on 
a given fixed renormalized mass is two-fold. In the static 
case it allows for a rapidly convergent iteration method, 
since the influence of cuts on the Green's function start 
only at p 2 = 9m 2 , where m is the physical mass. Sec- 
ondly, in the time evolution problem, the range of cutoffs 
which lead to results insensitive to the cutoff and the 
Landau pole value can be easily explored. Another re- 
cent development is the use of the 2-PI methods to obtain 
universal behavior at a critical point |ll| , where encour- 
aging results for critical indices have been obtained. 

The S-D equations we investigate here were studied 
earlier by Bender, Cooper and Guralnik 0, H3 and arc 
similar in structure to those obtained for the Gross-Neveu 
model and analyzed at all orders by Haymaker, Cooper 
et al. [2l]]. Our approach to the renormalization of the S- 
D equations parallels the treatment in that body of work 
and allows a simple renormalization scheme at order 1 /N 
where the renormalized vertex is replaced by 1. 

We organize this paper as follows. In Section [H] we in- 
troduce the auxiliary field formalism and discuss the 1 /N 
expansion as well as derive the unrenormalized S-D equa- 
tions for the two- and three-point functions. We also 
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discuss 2-PI expansion and the BVA. In Section II I II wc 
discuss the renormalization of the vacuum sector of the 
unbroken theory. In Section IIVI wc display the self- 
consistent renormalized S-D equations for the vacuum 
sector valid to order 1 /N. We solve these equations in the 
vacuum sector by an iteration scheme based on utilizing 
the lowest order in 1/N results and dispersion relations 
using a scheme used in Ref. Q. We compare our results 
to the related sunset approximation discussed in Ref. 3 
and find that at large coupling constant, g, there are 
significant differences between the sunset approximation 
and the next-to- leading order in 1/N self-consistent ap- 
proximation. Finally, in Section [V] we comment on the 
effect of the Landau Pole (triviality of continuum \<fi 4 
field theory j2<j) on our treatment of the 3+1 dimen- 
sional problem. 



A. The Large-N expansion 

The generating functional for the graphs of the auxil- 
iary field formalism is given by 



Z\j h K] = exp(iNW\j,K\) 



(2.6) 



= / dx JJd^cxpj i I dx C + ji(x)(j>i(x) 



+ iNK(x)x(x) 



The large-N expansion is obtained by integrating the 
Gaussian path integrals for 0j, letting each ji = j, 
and setting the free inverse propagator Go ij = <% Go 
(see H3,HE This results in an effective action 



II. AUXILIARY FIELD FORMULATION OF 
THE O(N) MODEL 

Consider the Lagrangian for O(N) symmetry: 



(2.1) 



S eff [x;j,K}/N. 



dx 



A 



2A 



\j°Goj+ l - TrlnGo 1 ^], 



K{x)x{x 

(2.7) 



where 

Go r [x]{x-y] = [n + x(x)]s(x-y), (2.8) 
and wc have introduced the notation 



Here, g denotes the scaled coupling constant g = X/N. 
The Einstein summation convention for repeated indices 
is implied throughout this paper. 

We introduce a composite field x( x ) by adding to the 
Lagrangian a term 



If— 9 ^ V 



2g { 2 L g 



(2.2) 



This gives a Lagrangian of the form 

fi 2 x(x) X 2 (x) 

9 2g 

which leads to the classical equations of motion 



(2.3) 



n + x (x) M*) = o-. 



and the constraint ("gap") equation 

x(x) = v + f 5>fro 



(2.4) 



(2.5) 



dx / dy ji(x) G[ X ]ij (x, y) j 3 (y) = N j o G o j . (2.9) 



The evaluation of the remaining path integral for x by 
steepest descent then leads to the 1/N expansion. The 
stationary phase-point of the integrand Xs[j, K] is deter- 
mined (implicitly) by the relation, 



K(x) 



X(x) + (J 2 - - j o G( , x)G(x, ) o j 



+ 2 G o(x,x) 



0. 



(2.10) 



Keeping only the Gaussian fluctuations we obtain for the 
W 



W\j,K] =W {0) + ^W {1) 



(2.11) 



=jXs ° ( y + y 2 ) + K ° Xs + ^ j o Go[xs} ° j 

+ l - Tr In Gq 1 [ Xs ] + ^ Tr In D^\j, K] 

where % s is to be viewed as a function of the sources j 
and K through Eq. 12.10f) above, and order 1/N 2 terms 
have been dropped. The bare inverse propagator for the 
auxiliary field x-i Do, is defined as the second derivative 
of the effective action with respect to x a t the stationary 
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phase point and its value is: 
1 



D- 1 \j,K](x,y) = --6 A (x,y) 



(2.12) 

N joG Q { ,x)G (x,y)G (y, )oj 



-G (x,y)G (y,x) 



The perturbative 1/N expansion for the connected 
Green's functions is obtained by treating all terms be- 
yond the Gaussian term in the effective action perturba- 
tively and is equivalent to a loop expansion in x for the 
effective action. Unfortunately this expansion has the 
same defect as ordinary perturbation theory when ap- 
plied to time evolution problems in that it displays secu- 
lar behavior as demonstrated in |15| beyond the leading 
order in large-N. It is precisely for this reason that the 
2-PI approach has proven so useful, since it leads to self- 
consistent S-D equation approximations that seem to be 
free from secularity. 

For completeness, we note that the generating func- 
tional for the 1-PI graphs, which is usually called the 
effective action, is the Legendre transform of W[j, K] to 
the new variables which are the expectation value of the 
fields 4> = SW/Sj, and X = 5W/5K. That is 

T[<f>, X }/N = W-jocl>-Ko X , (2.13) 

To order 1/N one obtains 

m x]/N = S cl [0, X ] + ^Tr In G^ 1 [ X ] + ^Tr In [<f> X ] , 



where 



(2.14) 



D- 1 ^ X ](x, y)=-~ $ 4 (x, y)-N cf>(x)G [x](x, y)4>(y) 
iN 

+ —G [ X ]{x,y)G[ X ]{y,x). (2.15) 



B. The S-D equations 

The S-D equations in the auxiliary field formulation 
treats the fields <f> and X and their two-point correlation 
functions on an equal basis. Thus for considering the S-D 
equations or the 2-PI effective action, it is useful to use 
the extended field and extended current notations, 

4>a(x) = [x(x),^i(x),^ 2 (x), . . . ,0jv(x)] , (2.16) 
j a (x) = [J(x),j 1 (x),j 2 (x), ■ ■ -,jN(x)} , a = 0,1, ...,N , 

where J(x) — NK(x), the source of field x introduced 
in the context of the large-N expansion. The generating 
functional Z\j\ and connected Green's function generator 
W[j] is given by the path integral: 



^ n <J 



,iS[<f>;j] 



(2.17) 



The path integral needs to be supplemented by stating 
the boundary conditions on the Green's functions. For 
initial- value problems one needs the CTP boundary con- 
ditions where in the time integration the time con- 
tour is a closed time path contour. However in discussing 
renormalization we only need to consider the vacuum 
equations and the Feynman boundary conditions on the 
fields. This is achieved by the usual ie prescription in 
deforming time slightly into the Euclidean region. The 
action S[<p;j] is given by: 

S[<t>lj] = - ~ f dx f dx (j) a (x) A~p[(j)}(x,x r ) (j)p{x) 



(2.18) 



+ / dx <f> a (x) j a (x) . 

Here, we have introduced the notation 





A-}[<j>](x,x')=( D ° 1( *> X,) 



(*0ij( X > x ') 



(2.19) 



with 



Do\x,x') 



5(x, x') 



GqIAx,x')= n + x (x) 5ij5(x,x') 



(2.20) 



The above are the diagonal entries in the Green's function 
matrix Gjj~*o[<£](x, x') defined as follows: 



(2.21) 



G o^M-- S(t)a (x)8Mz') 

= (d^{x,x') R^hxyy 

\K^{x,x') Gol&x'l 
while the off-diagonal elements are 
Kol[<j>}(x,x') = K^]{x,x') = ^(x)5(x,x'). (2.22) 
The S-D equations are obtained from the identity 

N r x 

n/ d ^^) elS[0uI = °- ^ 



0=0- 



The Heisenberg equations of motion and constraint de- 
scribing the time evolution of the O(N) model are ob- 
tained as 



E. 

fj 



(2.24) 

□ + x(x) <j>i(x)+Ki(x,x)/i=ji(x). (2.25) 
where the Green's functions G a p[j](x, x') arc defined by: 
54> a {x) S 2 W[j] 



G a/ 3[j](x,x') 



Sjpix') 8j a (x)5jp(x') 

D(x,x') K 3 (x,x'y 
Ki(x,x') G l j{x,x') 



(2.26) 
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Next, we introduce the 1-PI generating functional or ef- 
fective action T[cj>] by performing the Legendre transfor- 
mation 

= W\j\ - J dx<f> a (x) Ja (x) . (2.27) 
We obtain the equations of motion and constraint 



Sx(x) 



5<t>i(x) 



J{x) (2.28) 



9 

ji(x) 

D + x(x) (j>i{x) + Ki(x,x)/i 



(2.29) 



(2.30) 



We also define the inverse Green's functions 

D-\x,x') K~\x,x')\ 
y K-\x,x>) G^(x,x')) ' 

such that 

Ax"G-}[cj>}{xy)G M [i}{x",x') = 8 ay 8(x,x') . 

(2.31) 

The Green's functions G a p are obtained by inverting the 
equation 

G-}{x, x 1 ) = Golfa, x') + £ Q/3 (x, x') , (2.32) 

where Gq^(x,x') is given by Eq. (|2.21|) . and we have 
introduced the generalized self-energy matrix Yi a p{x,x') 
as 



33) 



By definition, the elements of the self-energy matrix are 
given as 



^00 


- Tl(x 


x') = 


1 5Gjj(x',x') 
2% S X (x) ' 


(2.34) 




-> Q,i(x 


x') = 


1 8G n (x',x') 
2i 54>i{x) 




S j - 


-> Qj (x 


x ') = 


1 SKj (x', x') 
i 6x(x) 




Sy — 


' ^ij 


x') = 


1 6Kj(x',x') 
i 8(j>i{x) 





The elements of the self-energy matrix are obtained by 
taking the functional derivative of Eq. H2.31f) . We have 



8G al3 [j}{xi,x 2 ) 
8<f>Jxs) 



- J dx 4 J dx 5 G a s[j](xi,X4) (2.35) 

X T S e 1 W\{x4,,X 5 ,X-i)G c p\j](x 5 ,X2) ■ 



Here, T a py denotes the three-point vertex function 

5G-i[<f>](x u x 2 ) 
T a/3j [(l)}(x 1 ,X2,x 3 ) = fj j2 - x (2.36) 



<50 7 (x 3 ) 



<5 3 r[^] 



Ha{Xl) 8(j)(j{x 2 ) 8^(x 3 ) 



In general, the three-point vertex function is the solu- 
tion of the exact equation 

5T, al3 (x,x') 



_ tL T a , 3l (x,x ,x ) = f a j3 7 8(x,x )S(x,x ) + 

J 
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<5</> 7 (x") 



(2.37) 

where f 0iiij = f ii0 j = f itjt0 = <% . The self-energy matrix 
depends implicitly on the three-point function as shown 
in Eq. i|2.35[) . From this equation, one can write a S-D 
equation for T in terms of an irreducible scattering kernel. 
We will derive this below in the unbroken symmetry case. 



C. Unbroken Symmetry Case 

These equations simplify considerably for the symmet- 
ric case, i.e. when 4>(x) = 0. The Green's function ma- 
trix is diagonal, Gu(x, x') = G(x, x') Sij. The self-energy 
matrix is also diagonal, Q,i(x,x') = £li(x,x') = 0, and 
Tiij{x, x 1 ) = S(x, x') 5ij. We have the gap equation 

X {x) = -H 2 + ^ G{x,x), (2.38) 
2,1 

and the Green's functions: 

D^ix.x 1 ) = - — S(x,x') + n(x,x'), (2.39) 
A 

G~ 1 (x,x') = \D + x(x)] 8(x,x') + E(x,x'), (2.40) 
where, by definition, the polarization and self-energy are 

II(x, x') = ^— J Ax\ J dx 2 G(x, xi)r(xi, x 2 , x")G{x 2l x') , 

(2.41) 

T,(x,x')=i J dxij dx 2 D(x, xi)T(xi, x 2 , x")G(x 2 , x') . 

(2.42) 

We notice that £ is of order 1/N since D is of order 1/N. 

In the symmetric case, there are S-D equations for the 
X&d vertex T which needs renormalization: Functionally 
differentiating the 4> inverse propagator with respect to 
x{ z ) w e can write 

V{x,y,z) = 5{x-y)8(x-z) (2.43) 
-iJdx l d X2 ^G{x-x 2 )D { x-x x )M { x u x 2 ,x^). 

where M is the 4> ~ X 1-PI scattering amplitude in the s 
channel (x\, x 2 ). In this paper we will use the schematic 
form 



1 - i DAIG . 



(2.44) 
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as shorthand for the above equation in either coordinate 
or momentum space when appropriate. 
The three graphs contributing to this are: 

M(x 1 ,x 2 ,X3,xa) = M stu (xi,X2,x 3 ,Xi) (2-45) 

+ J dx 5 dx 6 T(x 2 , x 3 , x 5 )D(x 5 ,x 6 )A(x 6 ,x 1 ,x A ) 

+ J dx 5 dx 6 T(xi,Xi, x 5 )G(x 5 ,x 6 )T(x 5 , x 2 ,x 3 ) , 

where 

. / \ SD~ 1 (xi, x 2 ) ,„ 
A( Xl ,x 2 ,x 3 ) = — ) u 27 (2.46) 

is the 1-PI 3-x vertex function which is finite, and 

M stu (x 1 ,x 2 ,x 3 ,x i )= 6T{X x 1 \ X2 \ X3) (2.47) 

M stu is the 1-PI in the s, t, and u channels scattering am- 
plitude for 4> — x elastic scattering. For renormalization 
purposes it is useful to have another representation of T 
in terms of K the 2 — PI in the s channel scattering ker- 
nel since this will facilitate the renormalization program 
needed below. 



D. BVA 

In the bare- vertex approximation (BVA) Q , the three- 
point vertex function T Q ^ 7 that appears in the S-D equa- 
tions for the inverse Green's function is approximated by 
keeping only the contact term, i.e. 



r a07 (x,x',x") = f a p 7 5(x,x')S(x,x") 



(2.48) 



This is justified at largc-N where the vertex corrections 
to the inverse Green's function equations first appear at 
order 1/N 2 . 

In this approximation, we have 



5G^(xi,x 2 ) 



1 (BVA) 



d(j)~ f (x) 



= - Gf ia {x\,x)f a p 1 Gp v {x,X2) , 

(2.49) 

which leads to the self-energies £ B J A given by 



^00 



II ^ (x, x ) — ^ G mn (x, x ) G mn (^, x ) , 



(2.50) 



(BVA) 



(BVA) 



(x , X ) — I Gi m (x , X ) Kjn (^1 x ) : 
(x, X ) = i K m (x , X ) G m j (x , X ) 



t(BVA) 



(x,x') 



Gij(x, x ) D{x, x) 
+ Ki(x,x') Kj(x,x') 



where we have used the symmetry property, Gij(x, x') = 
Gji(x',x) and Ki(x,x') = Ki(x',x). It can be verified 



by direct calculation, that indeed Cli(x,x') = fli(x',x), 
as expected. 

To summarize, in the BVA, one solves the equations of 
motion for 4>i{ x ) 

□ + X{x)] M x ) + K i( x , x )l l = > ( 2 - 51 ) 
and the gap equation for x{ x ) 

X (x) = -fi 2 + I {ti (*) + Gu {x, x)/i] , (2.52) 

i 

self-consistently with the equations for the Green's func- 
tions 

G-^x,x') = G-^x') + ^ A) (x,x') . (2.53) 

In the symmetric case the BVA polarization and self- 
energy are simply 



1 N 

rt BVA \x,x') = -^-G(x,x')G(x,x') 
£( BVA ) (x, x') = 1 G(x, x') D(x, x') . 



E. 2-PI expansion 



(2.54) 
(2.55) 



Here we review the 2-PI generating functional Q . We 
then derive the S-D equations that follow when we in- 
clude in T 2 the first term in a symmetrical (in propaga- 
tors x4>) 1 00 P expansion of the generator T 2 of the 2-PI 
graphs, namely the two-loop graph. To obtain a 1/N 
reexpansion, one then recognizes that in this loop ex- 
pansion, the three types of propagators have different 
dependence on 1/N (G ~ 0(1) and D,K ~ G(l/N)). 
The effective action is the twice-Legendre transformed 
generating functional: 



T[faG] = Sdas^H-TrOnlG- 1 ]} 

+ ^Tr{G o - 1 [0]G-l} + r 2 [G] 



(2.56) 



where 



6 2 S cl [fa 



6<j) a (x) 54>/s(x') 



(2.57) 



For initial- value problems, <f> is also a matrix in CTP 
space. S'dass is the classical Lagrangian written in terms 
of both <fi and x- I n the auxiliary field formalism the 
propagators for fab, <px and XX are treated on the same 
footing. Thus if we make a loop expansion the lowest 
order term in T 2 has two loops. In contrast, as discussed 
before, the 1/N approximation is asymmetric in x arL d ( t ) ^ 
since it is all order in <fi (for fixed X and a loop expansion 
only in X - 

The equations of motion for the field expectation val- 
ues follow by variation of T with respect to fa and x- We 
find 



ra 



X (x)] fa(x) = Ki(x,x), 



(2.58) 
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and 



(2.59) 



The equation for the two-point function follows by vari- 
ation with respect to G, which gives 



where 



2/ 



5G a p 



At the two-loop level we have that 



T 2 [G] 



1 

12 



Tr{ fGGGf} 

-\ [G lJ Gi J D + 2K l K J G l3 ] 



(2.60) 



(2.61) 



(2.62) 



Taking the derivatives of T 2 with respect to G a p, the £ 
matrix given by l|2.61[) leads to the same form for the 
matrix elements as we obtained previously from the S-D 
equations in the BVA. We notice, since D ~ 1/N and 
K ~ 1/N, the two contributions to T 2 are of order 1 
and 1/N, respectively, whereas the classical action is of 
order N. Thus if one is interested in counting powers of 
1 /N rather than loops, one would in next-to-leading order 
in 1/N ignore the second contribution to T 2 . This fact 
distinguishes counting loops in the auxiliary field formal- 
ism from counting powers of 1/N as noted in Refs. 0,0- 
The quantity ^[G] has a simple graphical interpreta- 
tion in terms of all the 2-PI vacuum graphs using vertices 
from the interaction term — \x$i4>i- When (4>) = 0, one 
obtains 



(2.63) 



^{m^-MH^ln^ 1 ]} 



+ - Tr^ 1 D + G„ 1 G 4> -2} + T 2 [G , D] . 

where G = {G^, D} and 

T 2 [G^,D] = (2.64) 
~\J d^i J dx 2 D(x 1 ,x 2 )G,p(xi,X2)G l p(x 2 ,xi) . 

The resulting equations for the two-point functions have 
no vertex corrections. 



III. ALL ORDERS RENORMALIZATION IN 
THE VACUUM SECTOR 

For renormalization it is necessary to only study the 
theory with unbroken symmetry since the renormaliza- 
tion is not changed when 0^0 (see e.g. Ref. jlcj)- The 



advantage of studying theory in terms of the auxil- 
iary field x ' 1S that the 1/N resummation improves the 
renormalizability as first discussed by Gross j^K ■ The 
renormalized theory can then be determined in terms of 
two "physical" parameters, which wc will choose to be the 
value of the (f> mass in the unbroken vacuum and value 
of the coupling constant at q 2 = 0. What we will find 
is that only the 4> two-point function needs wave func- 
tion renormalization and there is a Ward-like identity 
relating this renormalization constant to the vertex 
renormalization constant. The important effect of reex- 
pressing 4 field theory in terms of the x propagator D 
is that once the above renormalizations are performed, 
there are no further divergences in the elastic scattering 
of two (f> particles. 



A. Momentum space representation 

Here we consider the homogeneous case, relevant to 
understanding the vacuum sector. In Minkowski space, 
we define Fourier transforms as: 

G(x, x') = J [dp] e~ ip < x - x "> G(p) , (3.1) 
D(x, x') = J [dp] e-w-i*-*') D{p) , (3.2) 
where [dp] = > 2 J! d . Then, the gap equation is written as 

X = -M 2 + | J[dp}G(p), (3.3) 

and the equations for the Green's functions as 

N 



D-\p) 



X 



n(p), 



(3-4) 

G-» = [-p 2 + x] + E(p). (3.5) 
Finally, the polarization and self-energy are given by 
iN { 

H(p) = — J [dq] G(q) T{q,p- q) G(p - q) , (3.6) 
E(p) = l [ [dq] D{q) T(q,p- q) G(p - q) . (3.7) 



The only primitive vertex function for this theory is 
the vertex which we will just call T = Tuq. Keeping 
in mind the translational invariance properties, we write 

r(x,x',x") = J [dp}e-' lp < x - x "^ J [dq]e- iq < x '- x "^V{p, q) , 

(3.8) 

and the vertex equation 12.37f) becomes 

T{p,q) = 1 + AT(p,q). (3.9) 

The original S-D equation for T can be rewritten in terms 
of a 2-PI scattering kernel K 2 . which we schematically 
write in matrix form as 

T(p,p + q) = 1 + i[TDK 2 G](p,p+ q) . (3.10) 
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The primitive divergences of this theory have been dis- 
cussed in detail in Ref . [28| . The minimal degree of diver- 
gence (ignoring log improvements) is given by the simple 
formula 



D 



2B 



M. 



(3.11) 



where B is the number of external auxiliary field propa- 
gators D, and M is the number of external meson prop- 
agator (G) lines. Thus the meson propagator is naively 
divergent as A 2 and needs two subtractions (mass and 
wave function renormalization) , the D propagator is log 
divergent and needs one subtraction (coupling constant 
renormalization), and the vertex function has D= and 
is log log divergent. Another potentially divergent graph 
is the graph having four external meson lines, which in 
principal could be log divergent. However, the D prop- 
agators actually go as l/\n(p 2 ) at high momentum, and 
so the box graph with two D and two G propagators 
actually converges since 



d.r 



x(\nx) 



1 

InA 



(3.12) 



In this paper all the renormalizations will take place 
on mass shell, and we follow the treatment in |24|, since 
that approach is easiest to implement in the time evolu- 
tion problem, with the power series in — q 2 + m 2 becom- 
ing □ + m 2 acting on lattice versions of 5 functions. Our 
renormalization procedure will involve two steps. First 
we will identify the wave function and vertex renormal- 
ization constants as well as the physical mass of the 4* 
particle. We will then obtain naively finite equations for 
the multiplicatively renormalized propagators and ver- 
tex functions. Secondly we will show that when we re- 
place the bare vertices by the full vertex function minus a 
correction, we obtain finite renormalization constant free 
equations for the renormalized Green's functions and ver- 
tex functions. 

Before proceeding, let us define the multiplicative 
renormalization constants. We first identify the physi- 
cal mass by the zero of the inverse propagator for the <j) 
field, 



G 1 (p 2 = to 2 ) = 



(3.13) 



Once this mass is identified, then the wave function 
renormalization constant for d> is defined as 



Z^(m 2 ) 



dG-V) 



dp 2 



(3.14) 



If we write G 1 = — p 2 + x + E, then we also can write 
Eq. fl3~H|) as 



V(m 2 ) = 1 - 



dS(p 2 ) 



dp 2 



(3.15) 



The renormalized (f> propagator is then defined by 

G R (p 2 ) = Z? 1 G(p 2 ) . (3.16) 



The vertex function renormalization constant is de- 
fined by 



^r(PiP + q) = ZiT(p,p + q) 



(3.17) 



with the condition that on mass shell, with no momentum 
transfer. We have 



rfl(p,p)| p 2 =m 2 



which leads to 



Z7 l 



r(p,jj)|p2 



as( P 2 ) 



(3.18) 



(3.19) 



We will prove later that a Ward-like identity leads to 
the renormalization constants being equal, i.e. Z\ = Z-x, 
which tells us that the product TG is renormalization 
scheme invariant. 



B. Analyzing Divergences 

First let us realize that with our way of writing the 
Lagrangian, D(q 2 ) is a renormalization group (RG) in- 
variant and is just (apart from a sign) the running renor- 
malized coupling constant g r {q 2 ) for scalar meson scat- 
tering via single \ meson exchange. Thus g r (0) = g r is 
related to the scattering at q 2 = 0. In weak coupling, 
or in leading order in large-N, g r is the actual value of 
the scattering amplitude. However, in the full theory the 
scattering amplitude gets corrections from all the loops 
involving exchanging more and more D propagators. We 
have 



D-\q 2 =0) 
so that 



1 



1 



n(g 2 = 0) , (3.20) 



- = -- - ll(q 2 = 0) 

9 9r 



(3.21) 



If we wish to define a coupling constant renormalization 
constant via 
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Z g 1 g R 



then 



Z g = 1 + <7 R II(0) 



(3.22) 



(3.23) 



In terms of gn, the expression for the inverse propagator 
is 

D^iq 2 ) = -— +n[ subl V) , (3.24) 
9r 



with 



ni subl V) = n( g 2 )-n(o) 



(3.25) 
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Since n(q 2 ) is naively log divergent, then H^ suhl \q 2 ) is 
naively finite. At this stage we have that II = GTG, 
which is not yet in a form which displays its independence 
of the renormalization constants. 

Now let us look at the divergences of the <j> propagator. 
If we write G' 1 = -p 2 + \ + E, with E = GTD, and 



(3.26) 



then the bubble G(x, x) has quadratic as well as log diver- 
gences in 3+1, which are related to the mass and coupling 
constant renormalizations, respectively. In turn, the self- 
energy E has quadratic and log log divergences related 
to mass and wave function renormalization. Identifying 
the physical mass yields the relationship 







X 



E(p 2 



(3.27) 



To identify the naively finite part of E we now expand E 
around p 2 = m 2 

E(p 2 ) ee E + Ex(p 2 -m 2 ) + El sub2 V), (3.28) 

where E = E(p 2 = m 2 ), Ei = ^| (p 2 =m2) , and 
£[ sub2 l(p 2 ) i s naively finite and vanishes quadratically at 
the physical mass. Identifying 



Z^(m 2 ) 



dG-V 



dp 2 



(3.29) 



we can write 



G- 1 = 



so that 



P 



Z 2 E [sub2 



Z 2 E [sub2] 



) , (3.30) 



(3.31) 



This equation is now naively finite but not transparently 
independent of the renormalization constants. We will 
have to symmetrize the Dyson equation for E with re- 
spect to having fully dressed vertices at both ends in 
order to do this. 

Next we study the divergence of the vertex function. 
The vertex function can be written as 

T(p,q) = r(p,p)| p2=m2 + r[ subl l(p,g), (3.32) 

where r^ uhl \p,q) = T(p,q) - T(p 7 p)\ p 2 =m 2 is the ver- 
tex function once subtracted on the mass shell. Naively, 
the first term is log log divergent, and so rl sub L l (p, q) is 
naively finite. The renormalized vertex function on the 
mass shell with no momentum transfer is defined to be 
one, giving the equation 



rj?(p,p)| p 2 =m . 



or 



^r 1 



Zi T{p,p)\ p 2 =r 



r(p,p)| P 2= 



(3.33) 



(3.34) 



From this we obtain 

T R (p,q) = Z 1 T{p,q) = 1 + Z 1 T^{p,q). (3.35) 

If we write the S-D equation for T in the form T = 1 + Ar 
then we can rewrite this as 



T R (p,q) = 1 + Z 1 AT^ hl \p,q) 



(3.36) 



where the superscript means once subtracted on the mass 
shell. Again we need to show that the right hand side, 
when symmetrized appropriately, is independent of all 
the renormalization constants. 



C. Obtaining finite equations for the renormalized 
Green's Functions 

To remove the dependence of the above equation on 
Z\ and Zi the key ingredient is the symmetrization of 
the S-D equations, which usually have one bare and one 
full vertex function in their definition. We also need to 
prove the Ward-like identity that Z\ = Z 2 , which is done 
as follows: We have already said that D is RG invariant, 
so that the vacuum vertex function satisfies (for constant 
held X ) 



t(p,p) 



5G- l {p 2 ) 
5x 

Zf 1 T R (p,p) 



z; 1 



ZG r \p 2 
Sx 



(3.37) 



However, since x is n °t renormalized, we also have 

r*(p,p) = ^¥^. (3.38) 



Sx 

Thus we find that in this theory we have the identity 



Zi 



(3.39) 



This can also be shown by analyzing graphs. We have 



z x 1 = r(p,p)| p3=m 2 ee i 



dE( P 2 ) 



dx 



Z^ 



1 



<9E(p 2 ) 



dp 2 



, (3.40) 



(3.41) 



By studying graphs contributing to Z\ and Z-i the dif- 
ference between these graphs are seen to be naively fi- 
nite p8[ . and since these renormalization constants di- 
verge as log log A one has in the continuum Z\ = Z 2 . 
From this identity we find that the quantity TG is renor- 
malization scheme invariant and equals T R G R . We also 
have, Z\T = Z 2 T = T^. 

The next step is to get an integral equation for T which 
is iterative in T. The original integral equation for V that 
one gets by functional differentiation of the equation for 
the inverse x propagator is schematically 



1 



iDMG. 



(3.42) 
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where M is the one 1-PI in stu 4>X scattering amplitude 
defined in Eq. H2.45f) . We want to replace this equation 
by one in which the r.h.s also has a T. 

First we invert Eq. (|3.42|) to obtain the formal identity 
(see Ref. fZjfr) 



1 = T (1 - iDMG)- 1 



(3.43) 



We can therefore introduce a T into the integral equation 
using this identity: 

T =1 + iT (1 -iDMG)- 1 DMG (3.44) 
= 1 - iT DK 2 G 

Solving for K 2 and reintroducing the coordinates, we 
have 



M{x\,X 2 ,X 3 ,Xi) = K 2 {x 1 ,X2,Xs,X i ) 



(3.45) 



- J dx 5 da; 6 da;7dx 8 M(xi,x 2 ,x 5 ,x e )D(x 5 ,X7) 

x G(x 6 ,x s )K 2 (x 7 ,x s ,x 3 ,X4) , 

so that K 2 is the 2-PI in the s channel irreducible kernel 
of the Bcthc-Salpetcr equation for the scattering. What 
is important for renormalization is that the combination 
GDK is RG invariant as can be seen by looking at all 
the skeleton terms contributing to K. For example, if the 
skeleton one-particle exchange in the t channel contribu- 
tion to K , TGr, is investigated, we get the combination 
DTGTG, which is obviously equal to DrT r G r TG r . 

The key to eliminating the dependence of the renor- 
malized Green's function equations obtained earlier on 
Z\ and Z 2 is the symmetrizing of the self-energy and po- 
larization graphs with respect to the full vertex function. 
The strategy for doing this in quantum electrodynamics 
(QED) is found in the text book of Bjorken and Drell J2£j , 
where it is also shown that this procedure eliminates the 
problem of overlapping divergences. 

We have shown we can write 

r = 1 - i TDK 2 G = 1 + AT. (3.46) 
This equation allows us to substitute the bare vertex by 



1 = T-iTDK 2 G 



AT 



(3.47) 



The right hand side now multiplicatively renormalizcs 
exactly as T does, by the above argument. Replacing 
bare vertices by the right hand side will exactly absorb 
the extra factors of Zi and Z 2 . Equation 1)3. 46[) is useful 
to symmetrize E with respect to T. However a different S- 
D equation for T is needed when we symmetrize II which 
is a loop made of two <f> propagators. To obtain this S-D 
equation we realize that we can write D~ x as 

D~\xi,x 2 ) = [□ + X (xi)]S( Xl - x 2 ) + E[G, D](x u x 2 ) 

(3.48) 



where 



ST 2 \G,D] 

Eg, J (xx, a 2 )=-2i ; L J 3.49 
oG(xi,x 2 ) 



is just a function of the exact G and D propagators. 
Thus, using the chain rule and taking the total functional 
derivative of D~ x with respect to x(x3), we obtain the 
equation 



T(xi,X2,x 3 ) = S(xi - x 3 )S(xi 



x 2 ) 



dxsdxidx^dxe 



x r(x 3 , x 5 , x 6 )G(x 5 , x 3 )M 2 (x 3 , Xi , x 1 ,x 2 )G(x i ,x 6 ) 

+ A 3 (x 3 ,x 5 ,x 6 )D(x 5 ,x 3 )K 2 (x 3 ,X4, xi,x 2 )D(x 4 ,x 6 ) 

(3.50) 

In the above we have that A3 = 5D~ 1 /5x is the three % 
1-PI vertex function. The two scattering kernels are that 
M 2 is the 2-PI scattering kernel for cfxp — ► 4>(j> and K 2 
here is now the t channel kernel for \X ^ 4>4>- Explicitly, 
we have 



M 2 (x 1 ,x 2 ,xs,x i ) = 

K 2 (xi,X 2 ,X 3 ,X4,) = 



SG-^xux*) 
SG{x 3 ,x 4 ) 

5G~ 1 (xi,x 2 ) 
5D{x 3 ,Xi) 



(3.51) 



We note that in QED, A3 would be the three-photon 
vertex which is zero by Furry 's theorem |32j |. Again a 
study of graphs shows that Ar is rendered independent of 
renormalization constants when multiplied by Z\ = Z 2 . 
First, let us look at the Dyson equation for the renormal- 
ized vertex function. We use 



T(p,q) = 1 + AT(p,q), 



(3.52) 



where AT(p,q) = i [TGDK 2 ](p, q). Subtracting once on 
the mass shell leads to 

T R (p,q) = Z 1 T{p,q) = [1 + Z.AT^^p, q) 

= l + AT R (p,q), (3.53) 

which is now explicitly finite and free from any depen- 
dence on the renormalization constants, since 



AT R (p,q) = [T R G R D R K 2R }^ hl \p,q) 



(3.54) 



which is finite and written in terms of a renormalized 
skeleton expansion for the scattering kernel K 2R . Simi- 
larly using the second S-D equation we obtain 

AT R (p,q) =[T R G R G R M 2R ]^(p,q) (3.55) 
+ [A 3 ^^^ 2fl ][ subl ]( P , (? ), 

since D R - 0(1/N) and AT R (p,q) - 0(1/N). 

The equation for the x propagator contains the naively 
finite once subtracted polarization nl sub ^ = [GrG] [subl] . 
The right hand side of this equation however is not RG 
invariant. By inserting the second S-D equation for T in 
the form 1 = T — AT, then we have 



n [subi] = [( r _ Ar)GGr] [subl] 

-> [(rK-Ar^jG^Gflrfl]^ 1 ] 



(3.56) 
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and we see that all the dependence on the renormaliza- 
tion constants disappears. 

Next let us look at the (j> propagator. We have previ- 
ously shown that 

G R l (p 2 ) = m 2 -p 2 + Z 2 ^ suh2 ^p 2 ), (3.57) 

where £[ sub2 l(p 2 ) is the twice subtracted GTD and so is 
naively finite. So it is easy to see that once we replace the 
bare vertex by the first S-D equation for T in the form 
1 = r — Ar in S we again have 

G R \q 2 ) = ™ 2 -p 2 + 4 ub V), (3-58) 



This can also be written as 

to 2 = Z~L ml + S(m 2 ) - S(0) . 

For the vertex renormalization we now have 

9E(p 2 ) 



1(0) = T (P>P)\p 2 =0 = 1 + 



dx 



■> 2 =o 



(3.66) 



, (3.67) 



and the Ward-like identity is now Zi(o) = ^2(0) ■ 



IV. RENORMALIZED S-D EQUATIONS UP TO 
ORDER 1/N 



with 

E £ub2] = ^ S [ sub2] (3 5Q) 

= [(T R -AT R )G R D R r R }^ , 

and the subscript referring to subtracting twice on mass 
shell. 



D. Renormalizing at q 2 







In this subsection we would like to relate the physical 
renormalization discussed above by conventional renor- 
malization at the unphysical value q 2 = which is con- 
venient for evaluating integrals. For the inverse <f> prop- 
agator we now expand S(g 2 ) around the point q 2 = 0. 
That is, we let 



S(P 2 ) 



where now £ 



+ S 1(0)P 2 + £[T b V), (3.60) 



-0(0) 



0(0) 



S(p = 0), E 1(0 ) = 3^|(p2 =0 ), 



1— m, and 



Sg Ub2 '(p 2 ) is again naively finite, but vanishes quadrat- 
ically at p 2 = 0. In this case we identify 



7 -\ 

^2(0) 



dG'V) 

n.2 



dp 2 



p 2 =0 



(3.61) 



7- 1 

^2(0) 



V+ml + Z, 



2(0) ^0 



,[sub2]/ 2 



(3.62) 



We can write 
G~V) = 
so that 

G R \p 2 ) = -p 2 + m 2 + Z 2{0) 4 sub2l (p 2 ). (3.63) 

The same argument as before allows us to symmetrize 
the dependence of £ on the vertex function and obtain: 



4f V) 



(3.64) 



The finite mass m\ is related to the physical mass which 
is the zero of the renormalized Green's function via 







-TO + TOn + 



(3.65) 



Now that we have an all orders renormalization pro- 
cedure, to renormalize these equations at order 1/N we 
realize that Ar is of order 1/N so it can be ignored in 
the integral equations for II and E. Thus, it is consistent 
with the 1/N approximation to set = 1. This now 
gives us the finite renormalized equations for the inverse 
Green's functions, exact up to order 1/N. We can now 
write the renormalized S-D equations in the unbroken 
symmetry case as follows: 

For the 4> inverse propagator we have (Gij = G5ij): 



G~ R \q 2 ) = m 2 -p 2 + E^'V)- (4-1) 



sub 2]/ 2\ 



with 



V R (x,y) = iG R (x,y)D R (x,y) , (4.2) 
so that in momentum space 

E fl .(g 2 ) = i J[d d p]G R {q-p)D R (p). (4.3) 

We have 

4 ub V) =^(p 2 ) - E*(/=to 2 ) 

- (p 2 ~m 2 ) £ 1J? (p 2 = m 2 ). (4.4) 

For the % inverse propagator we have (g R = Xr/N) 

DrH<?) = -9 R l + n^V), (4.5) 

where 



n£ ub V) = iW)-n fi (o), (4.6) 



and 



iN 

^R(x,y) = — G R (x,y)G R (y,x), (4.7) 
or in momentum space 

K R (q 2 ) = ^ J[d d p}G R (q-p)G R (p). (4.8) 

In order to solve these equations, one first specifies the 
value of the renormalized mass, m 2 and the renormalized 
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coupling constant, g R . Then, both G R (p 2 ) and D R (p 2 ), 
and therefore T^ R (p 2 ), are just functionals of the finite 

quantity E^ ub2 '(p 2 ). Then Eq. (|4.4|l is the self-consistent 
equation we need to solved for the E^ ub2 '(p 2 ). 

One way of solving this equation at largc-N is to re- 
alize that E^ ub2 ' is finite and of order 1/N so that one 

can start with the value of E^ ub 2 ' found in leading order 
in the large-N approximation and iterate until conver- 
gence is obtained. The strategy for doing this in the 
related approximation of the keeping the three-loop sun- 
set graph in T 2 (which is the first term in a coupling 
constant reexpansion of the composite field propagator 
D) is very clear lydiscussed in the papers of H. Van Hees 
and J. Knoll 0, y| and we will essentially repeat their 
strategy here. Once E[ sub2 l is obtained, then one can re- 
construct G R and D R , and finally the subtraction terms 
E.r(to 2 ) and Ei^(to 2 ) needed for the renormalization of 
the time dependent evolution equations, as well as the 
renormalization constants Z\ = Z 2 - If one instead renor- 
malizes at q 2 = 0, one first specifies the parameters to 2 , 
and g R , and then gets a self-consistent equation to solve 
for Eg^ b2 '(p 2 ). Again one reconstructs the propagators, 
the self-energy and vacuum polarization, and determines 
the physical mass from the parameter to 2 , using 



jb 21 



J 0R 



(P 2 



(4.9) 



The advantage of the mass shell renormalization is that 
one definitely has a clear separation between the pole 
contribution to the dispersion relation for the Green's 
function discussed below and the three-particle cut which 
starts at 9m 2 . This large mass for the cut then allows a 
rapidly converging iteration strategy for E starting with 
the leading order in large-N propagators. 



A. General strategy 

Assuming that indeed the Green's functions G R (p 2 ) 
and D R (p 2 ) vanish at infinity (which is true in 3+1 di- 
mensions, but will need to be modified in 2+1 dimensions 
for D R (p)), then we can write the spectral representa- 
tions for the renormalizcd Green's functions as 



Gr( P 2 ) 



Dr( P 2 ) 



1 



-p 2 + ^ h2 \p 2 )- 
d(m 2 ) ImGflK) 
7r m\ — p 2 — ie 



~9r 



n 



ib 1] 



(p 2 ) 



d(m 2 ) lmD R (ml) 
it to 2 , — p 2 — ie 



(4.10) 

(4.11) 
(4.12) 

(4.13) 



where 



E fl (p 2 ) 



[d d q]D R {p-q)G R (q) 
° ^illm^(TO 2 ) 



(4.14) 
(4.15) 



n fl ( P 2 ) 



^l m G R ( m 2 )K^(p 2 ; m 2 , m 2 ), 

(4.16) 

ImG fl (m 2 ) (4.17) 



[d d q]G R (p~ q)G R (q) 

00 d(TO 2 ^ 



lmG R (m 2 )K^(p 2 ;m 2 ,m 2 ), 



with 



X[<V;to 2 ,to 2 ) 



[d d q] 



(m 2 — (p — q) 2 — ie){irv^ — q 2 — ie) 

[dV] 

(to 2 + {p E - <?i?) 2 )(m 2 + q%) ' 



(4.18) 
(4.19) 



with p 2 — > — p 2 E . The last (Euclidean 
for example in Itzykson and Zuber 
Feynman parametrization, as 



KW(p 2 



r(2-d/2) 
(47r) d / 2 



integral is found 
in terms of the 



(4.20) 



da [a(l 



a)p 2 + am 2 + (1 — a)m\ 



The advantage of writing things this way is now clear, 
as the kernel can be evaluated in d dimensions exactly 
The finite part of E^ can be obtained by performing the 
subtractions directly on the kernel at the physical mass. 
In 2+1 dimensions T, R requires one subtraction corre- 
sponding to mass renormalization, and in 3+1 there arc 
two subtractions related to both mass and wave function 
renormalization. II is finite in 2+1 dimensions and goes 
to zero at large momentum, which will lead to the need 
for using a subtracted dispersion relation for D R in three 
dimensions. In 3+1 dimensions II is log divergent so 
that one needs to subtract II once, which corresponds to 
coupling constant renormalization. These subtractions 
when done on the kernel K then automatically lead to a 
finite analytic (in the cut plane) expression for the sub- 
tracted kernel, which is then used to determine Et sub2 l 
and nl sub : 1 via the spectral representation. 

The above equations will be solved iteratively: We 
start with the free renormalized </> Feynman propagator 



G (s=p 2 ) 



1 



to^ — s — ie 
V 



(4.21) 



+ in S(s — to 2 ) . 
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FIG. 1: (Color online) AT| 4 r ' en (s; mf , m^) with mi = 
m — 1. 



FIG. 2: (Color online) ^2rcn( s > m i> m i) with mi = m.2 = 
m = 1. 



(For convenience, we shall drop for now the subscript i? 
in denoting renormalized quantities, and we will ignore 
the required subtractions for the polarization and self- 
energy.) Equation (|4.21|) . in conjunction with Eq. I|4.16|) . 

gives 



n (s) = \ f [d d q] G [(p-q) 2 } G (q 2 ) 

= - I&\s;m 2 ,m 2 ) . 
Using Eqs. I|4.11|l and (|4.14|) we obtain, in order, 
D (s) - 



(4.22) 



(4.23) 



So(a) 



^±1 Im£) ( m 2) KW{s;m\,m 2 ) . 

(4.24) 



We can now calculate the first update of the <f) propa- 
gator, G\(s), 

1 



Gi{s) 



m 2 — s + Eq(s) 



G (s) - AG(s) , 



with 



AG(s) = 



So(s) 



(m 2 — s) [m 2 — s + Eq(s)] 



(4.25) 



(4.26) 



The updates of the polarization and self-energy are ob- 
tained by combining Eq. (|4.25|l with Eqs. (|4.17l) and 
(|TJ5|) . We have: 

IL(s) =n (s)- r^HllmAG(m^l d l( S ;m?,m 2 ) 
Jo * 



1 [°° d(ml) 



2 Jo * 

d(m| 



ImAG(m^) 



(4.27) 



ImAG(m^)^ H (s;TO? 



2 ™2\ 



Ei(s) = E (s) 



d(m?) 



Im_Di(m^) 



(4.28) 



; r °° ImAG(m|)/v ^ (s; m?, m§) 



i] 



7T 



where 



Sib) = 



■<r 1 + n 1 ( P ) 



(4.29) 



We can now obtain the "new" correction AG(s), by re- 
placing Eo(s) with Ei(s) in Eq. (|4.26(l . Then, we iterate 
Eqs. H4.2fcip . I|4.27[l and (|4.28|l . until we achieve conver- 
gence. 



B. 3+1 dimensions 

The key ingredients for this case have been previously 
discussed in van Hees and Knoll . Those authors have 
used dimensional regularization arguments in order to 
regularize the kernel 



KW( S =p 2 ,mlm 2 ) 



d 4 q] 



(4.30) 



i\ — (p — q) 2 — ie](m\ — q 2 — ie) 



In the range (mi — m-i) 2 < s < (mi + 1TI2) 2 , one obtains 



K^(s,m 2 ,ml) 



16tt 5 



mi 



tan 



In 

-m 2 . 

m + m_ + s 



+ 7-2 + In 
A(s,m?,mf) 



mim 2 



47T/i 2 



A(s, mf , m 2 ,) 



— tan 



m + m_ — s 



O(e) 



A(s, m 2 , m 2 ) 

(4.31) 
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FIG. 3: (Color online) First iteration results in 3+1 dimensions (gn= 30). 



where 7 is the Euler-Mascheroni constant, and we have 
introduced the notation 

A(s, m\,m\) = J— [s - (mi + m 2 ) 2 ] [s - (mi - m 2 ) 2 ] . 

(4.32) 

In 3+1 dimensions, both the polarization, II(s), and 
self-energy, S(s), require rcnormalization. As advertised, 
the required prescriptions are 

n [subl](Q) = (4.33) 

S [sub2] (m2) =9sS [sub2] (s) | s=m2 = 0j (434) 

or 

ni subl ]( s ) = n(s) - n(o) , (4.35) 

s [ S ub 2](s) = s(s) _ s(m2) _ a sE(s) | s=m2 . (4. 36) 



As such, the renormalized polarization nl sub ^ (s) and 
self-energy £l sub2 l(s), will be calculated in terms of the 
kernels, 



K?l n (s; m 2 , m 2 ) = K^(s; m 2 , m 2 ) A^(0;m 2 , m 2 ) , 

(4.37) 



and 
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4 4 L( S ; m?, m 2 ) = KW(s; ml m 2 2 ) K^(m 2 ;m 2 , m 2 ) - (s - m 2 )d s K^{ S] mf, ml)\ _ 



We obtain 



^^(ssm^ml) - K[l cn {m 2 ;ml,ml) - (s - m 2 )d s K[> cn (s; m 2 , m\)\ 
[m\ — to 2 ,) 2 — s [m 2 + to 2 ) ( m\ 



I6n 2 s 



In 



A(s, mi, ?7i2 ) 



tan 



,s + m? — i s — m? + to 2 , 

' + tan" 1 



m 2 
l 



(4.38) 



(4.39) 



A(s, mf , m 2 . 



A(s, to 2 , to 2 ,) 



and 



167T 2 Af 4 S 



1 



167t 2 to 2 s 



167r 2 m 2 



m 2 A(s, mi, m2) 



tan 



m 2 

-™2 



(4.40) 



s + m? — m 2 



A(s, mi,m 2 ) 



tan 



s — m? + to 2 , 



A(s,mi,TO 2 ) 



S A(TO 2 , mi, 7712 ) 



tan 



m 2 + m? — 77i n 



A(to 2 , mi, to 2 ) 



tan 



_j m 2 — to 2 + 777 2 



A(to 2 , mi, TO2) 



(to 2 — to?,) 2 — m 2 (to 2 + m 2 ) 
A(to 2 , mi, TO2)to 2 



tan 



m + mf — m 2 

A(TO 2 , TOl, TO2) 



+ tan 



to — mf + to 



A(to 2 , TOl, TO2) 



Equations Q4.39|l and (|4.4C)|> arc analytically continued outside the range (mi — TO2) 2 < s < (mi + TO2) 2 , using 



tan 



_i s + 7n 2 — to 2 

A(s,TOl,TO 2 ) 



tan 



S — TO 2 + TO?, 



A(s,TOl,TO 2 ) 

1 ^ V (™l+"'2) 2 -s-^/ (mi— m 2 ) 2 — s 
1 \/ (mi+m 2 ) 2 — a+y/ (m 1 —m 2 ) 2 —s 



(4.41) 



In 



•J s— (mi— m 2 ) 2 — -^/s— (?ni+m 2 ) 2 



•y/ !-(mi-m!) 2 +\/ s— (mi+m 2 ) 2 



if s < (mi — TO2) 2 
, if s > (toi + TO2) 2 



We note that the final renormalized expression for K^] en 
could just as well have been obtained from covariant or 
non-covariant cutoff methods of regularization. 

For illustrative purposes, in Figs.nand[21 we depict the 
s-dependence of Kfj cn (s; to 2 , to 2 ) and K^j en (s; m 2 ,m 2 ). 
(We recall that Kf] en ( s; m 2 ,m 2 ) is related to the starting 
value of the polarization, Ho(s), see Eq. H4.22fl .~) 

We begin our iterations by calculating 



r 



Subsequently we iterate 
Gi(s) 



TO 2_ s + S [-b2] (s) 



G (s)-AG(s), (4.46) 



n 



[sub 1] 



r [sub 1] 
d(TO 2 ) 



( S ) = nr» 



(4.47) 



ImAG(TO 2 ) K [ *l n (s;m 2 ,m 2 ) 



<Jni- 



Go 00 



n 



[sub 1] 



1 



to^ — s — it 



K [ ^] en (s;m 2 ,m 2 ) 



D (s) 



9r 



n, 



[sub 1] 



,[sub 2] 



(s) 



d(TO 2 ) 



Im _Oq (TO 2 ) _/^ij^ cn (-^i TO 2 , TO 2 ) . 

(4.45) 



(4.42) 



(4.43) 



(4.44) 



d(TO 2 ) 



2 ./9m 2 n 



ImAG(TO 2 ) 



d(m 2 ) 



9m 2 



Di(s) 



ImAG(TO 2 )^ 4 r I cn ( S ;TO 2 ,TO 2 ), 
1 



][ sub2] (s) = ££ UUZJ (s) 



4m 2 1" 



+ nr i] ( S ) 

2]^ [°° d(mf) 

d(TO 2 ) 



(4.48) 



ImDi(m 2 ) (4.49) 



9m 2 



4m 2 T 

ImAG(m 2 ) K [ 2 4 l n (s;mlm 2 ) 
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The first iteration results are depicted in Fig. [3] Here, 
we compare our results with the leading order contribu- 
tion in the pcrturbative reexpansion of D(s) in powers of 
n(s), i.e. 



D(s) 



9R 



n(*) 
9 2 R [n 



(4.50) 



[sub 1] 



which is similar to the sunset approximation of Van Hees 
and Knoll. 

The effect of the self-consistent calculation is depicted 
in Fig.^J As stated by van Hees and Knoll Q| , the first it- 
eration results are negligibly modified by the subsequent 
iterations, since the main contributions come from the 
pole term of Gq , and the continuous corrections start at a 
threshold of s = 9m 2 . Henceforth, only the high momen- 
tum tail of the propagators is affected, and convergence 
is rapidly achieved. 
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FIG. 4: (Color online) Self-consistent values of the polariza- 
tion and self-energy in 3+1 dimensions. 
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FIG. 5: (Color online) (s; m\, m|) with mi = mi =m= 1. 




FIG. 6: (Color online) K^j en (s, m; m\, rri^) with mi = mi- 
m = l. 



C. 2+1 dimensions 

The 2+1 dimensions case is different from the 3+1 di- 
mensions case for two reasons. The first is that the vac- 
uum polarization goes to zero at large momentum so that 
D goes to a constant and obeys a once-subtracted disper- 
sion relation. Secondly, the self-energy requires only one 
subtraction to make it finite. However, we will find it 
more convenient for our iterative scheme to use two sub- 
tractions so that the renormalizcd Green's function has 
the same strength at the pole as the free one. 

Using Eq. I|4.2U[I . the kernel K^(s; m 2 , m|) is evalu- 
ated for the domain < s < (mi + m,2) 2 , and the result 
is analytically continued outside this range. We obtain 




FIG. 7: (Color online) First iteration results in 2+1 dimensions (<7=30). 



K^( S ;m 2 1: m 2 2 ) 



tan 



mi+m2 ' 
1 j n mi+m 2 + \/s 
87TVS mi +7712 — ' 

1 i n ■ v /s+(mi+m 2 ) 

8tt y/s [ (fni+m2) 



if s < , 

if < s < {mi + m 2 ) 2 
, if s > (mi + TO2) 2 . 



(4.51) 



In Fig.[S]we depict the s-dependence of i^ 3 l (s; m 2 , to 2 ), satisfy the boundary conditions required by the spectr 
which is related to the starting value of the polarization, 
n (s). We notice immediately that for large momentum, 
we have X[ 3 1(s;to 2 ,to 2 ) — > 0, like 1/y/s. In order to 
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representation, we use the simple manipulation 

D(s) = - g + [D(s)+g] (4.52) 

[°° d(m 2 ) ImJK) 
= - g + / 2 , 4.53) 

J 4m 2 it m A - s - te 

which gives 

£(s) = -ig I [d\]G{q 2 ) + ^ ^H)lmZ?(m?) 

J Jim 2 K 

^llmG(m 2 2 )K^( S -mlml). (4.54) 

7T 

Thus we see that one subtraction (mass renormalization) 
is sufficient to render the theory finite. However making 
only one subtraction, i.e. 



£[ subl l(m 2 ) 



. 



S [subi]( s ) = x(s) - £(m 2 ), 



(4.55) 



(4.56) 



one induces a finite wave function renormalization mak- 
ing the bare and renormalized Green's functions having 
different strengths at the pole. To avoid this, it is conve- 
nient to do a complete physical renormalization even in 
2+1 dimensions and instead consider: 



E [sub2]^ 



where 



d(ml) 



4m 2 

oo 



lm.D{m\) 



(4.57) 



d(m|). 



ImG^^^'cn^ m l. m 2) J 



K. 



L (s; to 2 , m^) = K^(s; m\, to 2 ) - K^(m 2 ; m\, to 2 .) 
- (s-to 2 ) d s K^{s-mlml)\ s=m 2 . (4.58) 

For the range < s < (mi + m 2 ) 2 , we obtain 



^2ren( S ; m l> TO 2) 



1 - 



1 

8^ 

111 



1 m 1 + m 2 + Js 
— In — 

sjs mi + m 2 - V s 

mi + to 2 + m 



2m 2 I m\ + m 2 — m 
(s — m 2 )(mi + rn 2 ) 



m 2 [(mi + TO2) 2 — m 2 



(4.59) 



The analytical continuation of the above result is done 
according to Eq. (|4.51(l . and we illustrate in Fig. [G] the 

s-dependence of K 2 3 J. en (s; m 2 ,m 2 ). 

For concreteness, we list the explicit equations we need 
the solve. For the first iteration, we have 



1 



Go («) = 



n„( S ) = il{l 3 l( S ;m 2 ,m 2 ), 



(4.60) 



(4.61) 



D (s) = 



1 



(4.62) 



4m 2 



(4.63) 



while the subsequent iterations provide the solution of 
the system of equations 

Gi(*)=- 1 [snh2] =G (s)-AG(s), (4.64) 



ni(«) = n (s) 

_ r d(m 2 ) 

J 9m 3 



(4.65) 



ImAG(m 2 ) if [3] (s ; to 2 , to 2 ) 



1 f 00 d(m 2 ) 



2 79m 2 1" 

<-~ d(m 2 ) 

lm 2 7 1 " 



Di(a) 



ImAG(m 2 ) 



ImAG(m^) K [3] { s; to 2 , to 2 .) 



1 



9- 1 + ni(a) 



(4.66) 



V,[sub2]/ s. ^[sub2]/ s. Z" 00 d(m 2 ) 2 

J (s)=XI J ( s )- / ImL>i(m 1 ) 

Jim 2 n 

<-°° d(m 2 



<hn 



2 7T 



(4.67) 

iy ImAG(m 2 )41j S ;m 2 ,m 2 ). 



Similarly to the 3+1 dimensions case, we plot the 
results after the first iteration (see Fig. UJ, and com- 
pare with the sunset-like approximation of van Hees and 
Knoll. Once again, we the corrections beyond the first 
iteration result are suppressed due to (mi +m 2 ) 2 thresh- 
old in the emergence of the kernels' imaginary parts, and 
self-consistent result virtually lies on top of the first iter- 
ation result. 



V. EFFECT OF LANDAU POLE 

What we have done earlier was a slight cheat for 3+1 
dimensions in that A</> 4 field theory is only an effective 
field theory in 3+1 dimensions, having nontrivial scatter- 
ing only when defined on the lattice (or with a momen- 
tum cutoff) , and the lattice spacing not taken to zero |29j | . 
The range of validity of the effective theory is determined 
by the position of the Landau Pole. The bare coupling 
constant A must be positive for the lattice field theory to 
be defined. Using Eq. (|3.21|) . we obtain the relationship 



fill 



fJ 



ffH(O) 



(5.1) 
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a = 



m 

5«n(o) ' 



(5.2) 



If we evaluate 11(0) in leading order, integrating over po 
and using a cutoff A in \p\ we have 



n(o) 



16tt 2 



p 2 dp 



(p 2 + TO 2 ) 2 



(5.3) 



The asymptotic behavior of the above expression at large 
cutoff, A, is 



ri(o) 



i 



16tt 2 



In- 



2A 



(5.4) 



In the cutoff theory one has that g^ is a monotonically 
increasing function of the bare coupling constant g, and 
has a maximum value defined by 



max 

9r 



1 



n(o) 



16tH 



l3 

TO 



(5.5) 



This behavior is shown in Fig. |SJ 

In order to capture all the physics of our approxima- 
tion, we would like the cutoff, A, to be much larger than 
the 9to 2 threshold, that is so important in getting the cor- 
rect physics. Thus, any A greater than say 30to 2 will be 
sufficiently large. From Fig. |H1 we conclude that as long 
as gn ~ 1; there is a wide range of momenta for which the 
effective theory is valid, and we can expect there exists 
a regime of cutoffs (less than the maximum momentum) 
for which the theory becomes cutoff independent. This 
behavior was shown to be correct in our mean field sim- 
ulations of disoriented chiral condensates j2^. In this 
regime, the continuum results we used here should offer 
a good approximation to the actual cutoff integrals re- 
quired for consistency in order that 4 field theory be 
a good effective field theory in the energy regime lower 
than the Landau cutoff for that coupling constant. To 
avoid Landau pole issues, using gn < 1 is much more 
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FIG. 8: The cutoff, A, dependence of the maximum value of 
the renormalized coupling constant, <?§' ax . 




FIG. 9: (Color online) Real part of the self-energy £o( s ), in 
3+1 dimensions, gn = 1. 
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FIG. 10: (Color online) Imaginary part of the self-energy 
Eq(s), in 3+1 dimensions, qr — 1. 



realistic than the value (gi?=30) chosen by van Hees and 
Knoll. In 3+1 dimensions at gn = 1, the resummed 1/N 
approximation and the sunset approximation are indeed 
very close. These results are illustrated in Figs. 151 andllUI 



VI. CONCLUSIONS 

In this paper we have discussed the rcnormalization of 
the S-D equations of the auxiliary field formulation of 
4> 4 field theory and then specialized to the self-consistent 
approximation to the coupled Green's function equations 
obtained by ignoring vertex corrections or equivalently 
expanded the 2-PI generating functional in loops and 
keeping the two-loop contribution to r 2 . We then ob- 
tained vacuum solutions for the self energy and vacuum 
polarization contribution to the inverse Green's functions 
in the bare vertex approximation. We compared our re- 
sults to the related sunset graph approximation of Van 
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Hees and Knoll and discovered that at strong coupling 
there were significant differences in these two approxi- 
mations. These differences become insignificant when gn 
is of order 1. In discussing the renormalization we con- 
fined ourselves to the case when there is no symmetry 
breakdown. When there is symmetry breakdown, one 
can use a "mass independent" renormalization scheme, 
such as that proposed by Kugo [34| , which requires only 
renormalizing the effective action in the symmetric the- 
ory, which is discussed here, to also renormalize the bro- 
ken symmetry case. Otherwise, one can study the more 
complicated set of broken symmetry Dyson equations and 
renormalize those as discussed in Ref. [2l|. The reason 
we chose physical on mass shell renormalization in this 
paper was that it allowed a clean iterative method for 
solving the self consistent S-D equations for the vacuum 
sector. For non-equilibrium uses of the effective action 
in the CTP formalism one docs not need to rely on self 
consistent solutions, since the time update equations will 
find these solutions dynamically. In the latter case it 
might be more efficient to determine the renormalization 
counter terms using the mass independent strategy of 
Kugo. 

An important question we have not addressed here is 
the fulfilment of Ward identities. In obtaining the renor- 
malized S-D equations, we relied on the Ward identity 
Z\ = Zi to obtain our final result for the exact renor- 
malized Green functions. After doing that we noticed 
it was consistent up to order 1/N to set the "renormal- 
ized" vertex function to 1. In leading order (and be- 
yond) in traditional large-N expansions of the auxiliary 
field generating functional, one explicitly can verify (as 



well as prove order by order) that Z\ = Z2, and also the 
Ward identities related to the O(N) symmetry, are satis- 
fied order by order. However, the resummation inherent 
in the self consistently determined 2-PI effective action 
can violate the exact Ward identities since one now has 
a result exact to a particular order in 1/iV, but differ- 
ent quantities (for example Z\ and Z2) may differ in the 
coefficients of the 1/N 2 and higher terms when the rele- 
vant Green functions are reexpanded in a series in 1/N . 
Therefore, although we have obtained finite renormalized 
S-D equations, which are exact to order 1/N one can- 
not rule out order (1/N 2 ) violations of Ward identities, 
Goldstone theorems etc. To obtain approximations in 
the 2-PI formalism which are consistent with Ward iden- 
tities and gauge invariance (for gauge theories) is still an 
unsolved problem and the subject of lots of recent atten- 
tion [3^, |3||. This subject is beyond the scope of this 
paper, which is mainly concerned with obtaining finite 
renormalized S-D equations for the truncated hierarchy 
of Green functions. In a future publication, we will use 
the results presented here to study the thermalization of 
renormalized 4 field theory in 2+1 and 3+1 dimensions 
for both N = 1 and the N = 4 O(N) model and study 
the grid sizes needed for the renormalized theory to be 
independent of grid size. 
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